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^ : ABSTRACT 

O 

We examine a collection of particles interacting with inverse-square two-body 

Q^ ! potentials in the thermodynamic limit. We find explicit large-amplitude density 

waves and soliton solutions for the motion of the system. Waves can be con- 

Q-i! structed as coherent states of either solitons or phonons. Therefore, either solitons 

>-Ch ■ or phonons can be considered as the fundamental excitations. The generic wave is 

- ^ , shown to correspond to a two-band state in the quantum description of the system, 

^^ 
j_j ' while the limiting cases of solitons and phonons correspond to particle and hole 

excitations. 
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1. Introduction and basic setup: There has been much recent interest in the 
Calogero-Moser-Sutherland model of interacting particles in one dimension [1,2,3] 
(which is often referred to in the physics literure as the CS model). This model is 
related to quantum spin chains with long range interactions between the spins [4], 
wave propagation in stratified fluids [5], random matrix theory [2,6] and fractional 
statistics [7]. 

The CS model is exactly solvable in both the classical and the quantum regime. 
Remarkably, the quantum solution is much easier to interpret, exhibiting a straight- 
forward analogy to the free fermion case. In a recent paper, Sutherland and Camp- 
bell examined the classical system in the thermodynamic limit and identified the 
excitations [8]. It was found that the classical system has solitons, corresponding 
to a single particle running through the rest of them, as well as small amplitude 
waves (phonons), identified with holes. The purpose of this paper is to derive 
large amplitude wave and soliton solutions of the classical system in the continu- 
ous limit, where the particles form a "fluid," and examine their correspondence to 
the quantum states. 

We consider a collection of particles with the hamiltonian 
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where for convenience we chose them of unit mass. In principle, such a system can 
be put in a box of length L (with an appropriate modification of the potential [2]). 
We shall be interested in the limit N,L ^ oo with N/L fixed. In this limit, the 
system can be described in terms of a density field p{x) and a velocity field v{x). 
At equilibrium, the particles will form a regular lattice of spacing a and density 
Po = l/a. The particle current is J = pv and by particle conservation 

p + dJ = p + d{pv) = (2) 



where d = d/dx. The kinetic energy of the system is 



K = I dx^pv 



We can formally solve (2) for v to obtain v = —d ^p/ p, and the expression for the 
kinetic energy becomes 

This is exactly the kinetic term of the collective field hamiltonian description of a 
many-body system [9]. The potential energy can also be expressed in terms of the 
density. The naive expression, however, which would be 

is incorrect. The reason is that the interaction is singular at coincidence points, 
and thus a substantial part of the potential energy comes from nearest neighbors 
and is not accurately reproduced by the naive continuous expression. The cor- 
rect expression requires a careful conversion of the discrete sum in terms of the 
continuous fields. Alternatively, we can simply take the classical limit {h -^ 0) of 
the quantum mechanical expression derived in the collective field formulation [10]. 
The result is 

where p stands for the Hilbert transform: 
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p= dy P.P. piy) = i hm / dy ( + -) piy) (5) 
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The first term, which accounts for the interaction of each particle with its few 
nearest neighbors, is the dominant one in the limit where the scale of variation of 
p is much larger than the lattice spacing. In our case, however, we are interested 
in finite-width fluctuations and we must keep the full expression. 



The dynamics of the system can be found by varying the lagrangian L = 
K — V + fip with respect to p. The chemical potential fi plays the role of a 
Lagrange multiplier ensuring that the total number of particles remains constant. 
The resulting equations of motion are 

fi-i- 12 ^^9 2 , o~ , 9fdp\^ 9r.(dp\ 
-d v-^v -—P +9 dp + -\^—j +-(9(^— j +/i = (6) 

as well as (2). The inverse derivative operator in (6) is defined in terms of the 
principal value, in Fourier space d~^ = lim£_+o k/{hP' + e^). In particular, acting on 
a constant it gives zero. By requiring that the static configuration v = 0, p = po 
be a solution of (6), we obtain the value of the chemical potential 

/^ = —Po (7) 

This is in agreement with the value obtained from the exact solution of the many- 
body problem [2,8]. 

2. Small- amplitude waves: From the above equations we can obtain the dis- 
persion relation in the linearized regime of small-amplitude waves, which we shall 
call phonons. Noting that the Fourier transform of dp is 7r\k\p{k), we obtain 



Vphase=(^-^j =9[^Po-y) °^ ^ = V9{^Po\k\- — 



From (8) we deduce that the velocity of sound Vg, defined as the phase (or group) 
velocity in the long wavelength limit, is 

Vs = ^^po^/9 (9) 

In terms of the group velocity Vg the dispersion relation becomes 

^ = -^^ (10) 

We observe that (9) and (10) are the exact results. The group velocity is always 
smaller that the velocity of sound, and the above linearized waves can be identified 



with holes in the quantum theory. Notice that the above formulae are valid for 
|A;| < vrpo = ^, else the group velocity turns negative. This is reasonable, since the 
above condition restricts the momentum to the fundamental region of the Brillouin 
zone, thus avoiding umklapp. 

3. Solitons: As observed in [8], the many-body system should exhibit soliton 
solutions, corresponding to particle excitations. On the other hand, in [11] an 
equation similar to (6) was written for a system of free fermions, coming from 
an effective lagrangian chosen so as to reproduce the full quantum mechanical 
dispersion relation of the system at the semiclassical level. This equation has 
solitary wave solutions [11]. As we will demonstrate here, our equations (6), (2) 
also have solitary wave solutions of a rational type; we shall call these solutions 
solitons, and will comment later on their true nature. 

For a localized constant profile configuration, propagating at speed v, both p 
and V are functions of x — vt only. From (2) we have 

divp — vp) = and thus v = v (11) 

P 

In the above, the integration constant is fixed by the boundary condition that 
w — * at a; ^ ±oo, where p -^ po. Similarly, (6) becomes 

To guess a solution for (12) of the form psoi = Po + Sp, where 6p is localized, we 
notice that the term in (12) containing the Hilbert transform will always produce 
out of a localized function a tail falling off quadratically. Thus, 6p itself should 
have such a behavior at infinity. The simplest function of this form is 



Psol = Po + 
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Plugging the above form in (12) we find, after an amount of algebra, that it is 



indeed a solution, provided that v > Vg and 
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We finally arrive at the soliton profile 
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(13) 



(14) 



The above solution is, strictly speaking, a solitary wave. True solitons are solitary 
wave solutions of integrable equations, and scatter off each other preserving their 
number and asymptotic momenta. Since the initian many-body system (1) is 
integrable, we expect the corresponding continuum system to be also integrable, 
although a direct prrof is lacking, and thus (14) to be a true soliton. This is 
corroborated by the correspondence of these solutions to particles, as demonstrated 
below. 

The above soliton carries particle number Q, momentum P and energy E, 
defined as the extra amount over the static solution po- We find 



Q = dx (psol - po) = 1 

P = dx Psol 1^ = V (15) 

E = jdx [ir(psol) + V(psol) - Vipo)] = iv2 



We observe that the net particle number carried by the soliton is 1, independently of 
its velocity; its momentum and energy are also those of a free particle of unit mass 
moving at the soliton velocity v. Therefore, the soliton can be exactly identified 
with a particle excitation of the system. Again, this is in agreement with exact 
results drawn from the quantum theory, where particle excitations always move 
faster than the sound [8]. Notice, further, that the solitons become thinner as their 
velocity increases, while their spread diverges as they slow down to the velocity of 
sound. 



The above result for Q implies that the displacement of the equilibrium lattice 
far away from the soliton is ± half lattice spacing either way (so that there is 
an excess of one particle near the soliton). This result, as well as the form of the 
soliton (14), are at odds with the results found in [8]. We suspect that the source of 
the discrepancy is the truncation to a finite number of a;-derivatives of the form for 
the potential in [8]; this turns the equation to a local one and gives the soliton an 
exponential decay, rather than the inverse-square decay of the nonlocal equation. 
We also notice that our soliton has some important qualitative differences from the 
solitons in the semiclassical fermion theory of [11]: Our solitons carry a positive 
particle number of 1, as opposed to a negative particle number in [11], which would 
rather identify them as holes. Further, there are no static solitons in our case, since 
|v| > Vs, while in [11] solitons can slow down to zero speed. Finally, the definition 
of momentum used in [11] differs from ours by a surface term. Clearly (15) is the 
physically sensible definition in our case. 

4- Finite amplitude waves: Soliton profiles moving at very large distances from 
each other will obviously remain solutions. If we could form a state consisting of 
a sequence of solitons at regular distances spaced by A, all moving with the same 
velocity v, we would have found a large-amplitude wave solution with wavelength 
A. We thus try the form 

^ / \ 1 sinh^ 

Pwave - Po = > Psol(a; - nX) - Po) = t- —^ 2',x (^6) 
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where now the parameter u is not necessarily given by v^/(v^ — v^), since the 
proximity of the other solitons may have changed their common velocity. The above 
waveform is characterized by its amplitude A, defined as midway the distance from 
peak to trough, 

^^Pmax^Pmin^ 1 ^^^^ 
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as well as by its wavelength A. Substituting the form (16) in (12) we find, again 



after quite a bit of algebra, that it is indeed a solution provided 
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The above is the amplitude-dependent dispersion relation for the nonlinear waves 
of the system. Before we interpret it, however, we must note the following: The 
conventions used for deriving (12) were that the solution p carries some particle 
number and momentum on top of the "vacuum" solution po- This is reasonable for 
an isolated soliton, but rather inconvenient for a wave solution, which is thought 
to be a fluctuation carrying no net particle number and no net momentum (no 
drift). But the presence of the solitons in (16) adds one particle per length A, and 
thus the true equilibrium density of the system is po + j- Further, the solitons 
contribute a momentum v per length A; to neutralize it, we must boost the whole 
system in the opposite direction by an appropriate amount. After performing these 
redefinitions, the expression for the wave in terms of the true velocity v and true 
background density po is 

A yVA2A2 + 1- AAcos^ J 
and the nonlinear dispersion relation in terms of the amplitude A is 
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In the limit A ^ oo the above equations reduce to the single soliton solution. 
In the limit A ^ 0, on the other hand, the above formulae become 

Pwave — Po ~'t~ A COS KX , K — —— 

r- (21) 

V = — = Vs — K 

k 2 

which is the small amplitude wave solution and dispersion relation. We see, there- 
fore, that the above solutions interpolate between the two extreme cases. We stress 



that the generic wave can run either faster or slower than the speed of sound. It 
should also be noted that the above wave solution constitutes a solitary wave for 
the continuum limit of the system with periodic space (that is, the Sutherland 
model), where the period is the wavelength. 

5. Discussion and conclusions: In summary, we have found exact soliton and 
wave solutions for the CS system in the continuum limit. Certainly the above 
do not exhaust the list of solutions; the general motion of the system will be a 
nonlinear superposition of waves (or solitons). Although we could find such many- 
soliton or many-wave solutions, it is an algebraically laborious task of not much 
interest. It serves, nevertheless, as an indication that the above solitary waves are 
true solitons, as expected from the integrability of the original model. 

It is instructive to put the above solutions into correspondence with the quan- 
tum mechanical states. Consider A^ particles in a space of length L. The ground 
state of the system consists of a "Luttinger sea" in the pseudomomentum, with 
spacing between adjacent particles equal to 1-nljL and "Fermi level" TidN/ L, where 
g = i{i — h). At the limit h -^ 0, N,L —^ oo, N/L -^ po, the ground state becomes 
a continuous filled band with Fermi level Pp = n^ygpo- A small amplitude wave, 
corresponding to a hole, is a very small gap in the band. A soliton, corresponding 
to a particle excitation, is a single particle peeled from the Fermi level and placed 
some distance above. The generic finite amplitude wave corresponds to a state 
with two continuous filled bands, of widths Pi and P2 (with Pi + P2 = 2Pp) and 
with a gap G between them. These are related to the wave parameters as 

^^ _ 2TVg 

pfV G ^ '22> 
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Such a state can be visualized as arising either by successively exciting single par- 
ticles by the same constant momentum, until they form a continuous band, or by 
gradually augmenting the gap of a hole, until it becomes finite. This state can 
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thus be thought of as either a coherent state of sohtons (much hke the way we 
constructed the wave solution), or as a coherent state of phonons, their nonhnear 
nature accounting for the change in profile as they accumulate. Indeed, the soli- 
ton itself can be thought of as a superposition of many phonons with very large 
wavenumber, and the phonon as a superposition of many solitons just above the 
Fermi level. For the finite N (finite L) system the distinction between the two is 
fuzzy and in principle only one kind of excitations need be considered as fundamen- 
tal. Note, further, that quantum mechanically the holes behave as particles with 
fractional statistics of order h/i (meaning that i/h of them put together would form 
a fermion). At the classical limit h ^ 0, thus, they become bosons, as they should 
be since phonons obey no exclusion principle. Particles, on the other hand, carry 
statistics of order i/h. Thus in the classical limit they become "superfermions" 
meaning that no two of them can occupy relatively nearby quantum states. This 
is consistent with the inverse square repulsion between the classical particles. 

The above results are of direct relevance to the large- A^ limit of one-dimensional 
free matrix models. The particular wave and soliton solutions correspond to mo- 
tions of the density of eigenvalues in the unitary and hermitian models, respectively. 
Taking, for clarity, the hermitian case, the motion of a free N x N matrix M with 
angular momentum i is 

Mjk = Sjk (Pjt + aj) + {l- 5jk)—^ (23) 

Vj Pk 

The situation where most of the eigenvalues lie on a regular lattice with only one 
of them moving with velocity v is reproduced by choosing 



Pj 



27r£ / . A^ 



a{N -2)V 2 



iV \ 
J - j) (for J < AT) , pN = y, aj = (24) 



(Notice that the above momenta Pi, ■ ■ -Pn-I span the values between the two 
"Fermi" levels ±^-) It should be possible to prove analytically that the eigenvalues 
of (23) with parameters (24) have a density as given by our soliton solution, but in 
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practice this is a very hard task. The corresponding problem for unitary matrices 
is even harder to tackle, while our wave (19) readily provides the solution. Many- 
soliton solutions will be given by eigenvalues of (23) with, now, more than one of 
the momenta pj taking values equal to the velocities of the solitons, while the rest 
span the "Luttinger sea". 

The solutions found in this paper are very similar to the ones in stratified fluids. 
This is sensible, since the motion of these fluids (under some conditions) is governed 
by the Benjamin-Ono equation, which is known to have solitons behaving like 
Calogero particles [5]. This also suggests that the the many-soliton solutions of the 
CS model will correspond to the ones of the Benjamin-Ono equation, at least when 
all of them are left- or right-movers. The interesting fact is that stratified fluids 
themselves behave, in this respect, as hydrodynamic collections of CS particles. 
The exact mathematical connection of the two systems is still obscure. 

We conclude by noting that the quantum mechanical problem separates into 
two noninteracting chiral sectors, having to do with excitations near either end 
of the Luttinger sea. (The two sectors mix nonperturbatively when a number 
of particles of order A^ is excited, depleting the sea.) Therefore, the equation 
(6) governing the continuum system should also decompose into two nonmixing, 
first-order in time equations, one for each sector. For the corresponding equation 
for free fermions this is indeed the case [12]. In fact, from the collective field 
description of the system when only one chiral sector is present, we deduce that 
the chiral equations are exactly of the Benjamin-Ono type [10,13]. The exact field 
combinations in terms of which this decomposition would be achieved, however, 
are not known and constitute an open problem. 
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